Ground State of a Heisenberg Chain with Next-Nearest-Neighbor Bond Alternation 
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We investigate the ground-state properties of the spin-half J1 — J2 Heisenberg chain with a next- 
nearest-neighbor spin-Peierls dimerization using conformal field theory and Lanczos exact diago- 
nalizations. In agreement with the results of a recent bosonization analysis by Sarkar and Sen, we 
find that for small frustration (J2/J1) the system is in a Luttinger spin-fluid phase, with gapless 
excitations, and a finite spin-wave velocity. In the regime of strong frustration the ground state 
is spontaneously dimerized and the bond alternation reduces the triplet gap, leading to a slight 
enhancement of the critical point separating the Luttinger phase from the gapped one. An accu- 
rate determination of the phase boundary is obtained numerically from the study of the excitation 
spectrum. 
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Asymmetric spin-half ladders have recently attracted 
much theoretical interesliiiSi^*^ due to possible ex- 
perimental realizations in delafossite cuprates such as 
YCUO2.5 (Ref. 0) and the unusual physical effects that 
the asymmetry in the leg exchanges could introduce both 
in the ground-state correlations and in the properties of 
the excitation spectrum. The simplest and most trans- 
parent Hamiltonian representation of these systems can 
be given in terms of frustrated {J1 — J2) spin-half Heisen- 
berg chains with a spin-Peierls dimerization in the next- 
nearest-neighbor (n.n.n.) interaction, namely 
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where Ji , J2 > are the nearest- and next-nearest neigh- 
bor exchange couplings, respectively and < 5 < 1 is 
the bond alternation parameter. In the following, all the 
energies are expressed in units of Ji . For (5 = the prop- 
erties of the model and the existence of a quantum phase 
transition at ( J2/ Ji)c — 0.241 are well establishediL^ Be- 
low the critical point the system is in a Luttinger spin- 
fluid phase with a gapless spectrum, no broken symme- 
try, and a finite spin-wave velocity. Spin correlations are 
characterized by a power-law decay defined by the same 
exponent of the pure Heisenberg chain (J2 = 0). In the 
regime of strong frustration, instead, the ground state is 
spontaneously dimerized, two-fold degenerate and adia- 
batically connected to the Majumdar-Gosh exact solu- 
tion for J2I Ji = 0.5Siifl with exponentially decaying spin 
correlations. 

The effect of the bond alternation on the n.n.n. 
interaction has been recently investigated using field- 
theoretical approachesi*^ with remarkable disagreement 
in the conclusions. In fact, Chen and collaborators^ us- 
ing bosonization and a renormalization group analysis 
claimed that in the limit of small frustration the asym- 
metry in the n.n.n. integrals destabilizes the isotropic 



Heisenberg fixed point leading to a new phase with gap- 
less excitations and vanishing spin-wave velocity. On the 
other hand, with a similar bosonization approach Sarkar 
and Sen^ found that bond alternation represents an ir- 
relevant perturbation in the regime of small frustration. 
In this paper, using conformal field theory we show that 
the n.n.n. spin-Peierls operator, i.e., the one associated 
with the alternation in the n.n.n. exchange, 
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represents indeed an irrelevant perturbation for the 
Heisenberg chain in the regime of weak frustration. Fur- 
thermore, using Lanczos diagonalization technique we 
find that in the dimerized phase the effect of the bond 
alternation is to reduce the spin gap leading therefore to 
a slight enhancement with 5 of the critical value { J2/ Ji)c 
for the fluid-dimer transition. Finally, we present an ac- 
curate numerical determination of the shift introduced 
by the bond alternation on the phase boundary between 
the Luttinger and the dimerized phase. A preliminary 
account of this work was already presented in Ref 

In general, the relevance of the operator ||2Jl for the 
Luttinger fixed point can be determined by calculating 
its scaling dimension X. This can be defined in terms of 
the size-dependence of its generalized susceptibility 

Xo = j{^o\d{Eo - n)-^d\^o) « L^d-^) , (3) 

where L is the number of sites of the ring, so that A" < 1 
characterizes a relevant operator pi2i Here \iPq) represents 
the unperturbed Luttinger ground state. The calculation 
of X can be done by identifying the quantum numbers 
{j} (referenced to the ground state with energy Eq) of 
the intermediate states |j) appearing in the Lehmann 
representation of xo ^ In fact, in a conformally invariant 
field theory, the scaling dimension A of a given operator 
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is related to the finite-size corrections of tfie energy of 
the lowest intermediate eigenstate \ j) by 

AE{L) = Ej{L) ~ E(i{L) = 2ttv sX/L , (4) 

where Vs is the spin-wave velocity. In the case of the 
n. n. n. operator ^ , the quantum numbers of the inter- 
mediate states, referenced to the ground state, are the 
total spin S = 0, momentum k — n and even parity un- 
der spatial reflection {I ^ L — I), R = 1. In all cases, 
the only difference with the nearest-neighbor (n.n.) spin- 
Peierls operator, 

On.n. = ^(-l)'Si • S, + i , (5) 

I 

the well-known relevant perturbation, is the spatial- 
reflection quantum number, R = —1. The finite-size 
corrections can be computed cither by Bethe ansatz (for 
J2I Ji = 0) or by bosonization. According to the stan- 
dard notation, at the Luttinger fixed point, the effective 
Hamiltonian describing spin fluctuations is given by a 
free boson thcory>i^ 

Hl = yJ {^n^(^) + ^[dM^)?] , (6) 

where Vs is the spin velocity, is the spin fluctuation 
field and n(a;) is the conjugated moment. The dimcn- 
sionless coupling constant K is the key parameter of the 
bosonization theory, whose value can be analytically de- 
termined at J2 = by comparing the low-lying excitation 
spectrum of the Hamiltonian © with the Bethe ansatz 
solution of the lattice model: K = {2 — 2cos~^ A/tt)"^. 
Here < A < 1 is the coupling constant defining the 
easy-plane anisotropy: at the spin-isotropic point A = 1, 
leading to K = 1/2. The finite-size corrections of the 
lowest excitation energies of this Hamiltonian arc given 
byOi 

£;„,™(i) - E,{L) = ^ + . (7) 

For each pair of integers n and m {n = m ~ being 
the ground state), a low-energy excited state with en- 
ergy En^m is defined by choosing appropriate boundary 
conditions for the fields: 

$(X)-$(0) = nV^, (8) 

L 

J dxU{x) = 2mV¥. (9) 


Using the relation between the boson fields and the phys- 
ical electron operator ones, it was then derived in Ref.flSl 
that n is the z component of the spin, while m -|- n is the 
momentum of the low-energy excitation in units of tt. 
Following this work, it can be shown that the reflection 
symmetry TZ acts on the fields as follows: 

7^t$(a;)7^ = -$(L- x) , (10) 
7^t^(a;)7^ = -n{L-x) , (11) 
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FIG. 1: Size scaling of the gap of the lowest eigenvalues of 
the isotropic spin-half J1 — J2 chain with 5*^ =0, k — n, even 
parity under spin reflections, and R — ±1 under lattice reflec- 
tions: J2/J1 ~ (circles), 0.1 (triangles), and 0.2 (crosses). 
Lines are guides for the eye. 

implying that the quantum number m changes sign under 
spatial refiection TZ, whereas n does not, and, therefore, 
the ground state has R = 1. For m ^ 0, each pair ±m 
leads to two degenerate eigcnstatcs with opposite reflec- 
tions, R ~ ±1. 

Clearly, the excitations of Eq. ^ do not exhaust the 
full energy spectrum of the Luttinger model because, on 
top of them, the usual boson spectrum is also present: 

Ek{L) ^ VsY,\k\aW , (12) 

where aj. is a boson field defined in terms of <I>(x) and 

H(a;), which satisfies iV' a\jZ ~ ^^-a;- 

For both the spin-Peierls operators previously defined, 
the appropriate choices of quantum numbers are n = 
(i.e., 5*^ = 0) and m = odd (i.e., momentum tt). The 
two lowest excitations with m = ±1 and n = can be 
combined to construct two states with defined reflection 
parity, R = ±1. In the spin-isotropic case, K = 1/2, the 
excitations with (m = ±1, ri = 0) and (m = 0, n = ±1) 
are degenerate, the latter ones corresponding to a i? = 1 
triplet (the spatial reflection does not change n). There- 
fore, the state with n = Q and i? = 1 is the 5^ = state 
of the triplet, and the remaining state with n = and 
R = —1 has to be a singlet. Analogously, all finite-m 
excitations with n = are forbidden to be singlets with 
i? = 1 by the spin isotropy. In order to obtain a low- 
energy singlet excitation with i? = 1, we have to com- 
bine the m 1 (i? = — 1) excitation with at least three a\ 
bosons with zero total momentum. In fact one boson can- 
not have zero momentum and two bosons cannot change 
the spatial reflection symmetry. The minimum energy of 
such state is readily evaluated [with fci = fc2 = 27r/L, 
and fcs = — (fci -I- fc2)] as 

A£:(i) = 2t:Vs{K + 4)/L , (13) 
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FIG. 2: Size scaling of the triplet gap at J2/J1 = 0.5 for 5 = 
(full circles), 0.25 (squares), 0.5 (triangles), 0.75 (crosses), and 
1.0 (empty circles). Lines are quadratic fits. Inset: thermo- 
dynamic gap versus S obtained by extrapolation (sta rs and 
dashed line) and with the variational approach of Refs. 
(continuous line). 

leading to the result X ~ K + 4 for any operator with 
the same symmetries of the operator 

Although the previous equation has been extracted in 
the isotropic case (A = 1), where K = 1/2, it is indeed 
valid in general for < A < 1, as can be verified, analyti- 
cally, in the XY model (i.e., at A = 0) or, numerically, for 
A < 1 (not shown). The reason is that, in the entire gap- 
less regime, the low-energy excitations are described by 
Eqs. O and ifT^ . Notice that, for A 7^ 1, the total spin is 
no longer a good quantum number, nevertheless the par- 
ity under spin-reflection [(5^, Sy, 5^) ^ (-5'^, , -S'^)] 
can still be used to label the low-lying states. 

The previous analysis shows that the finite-size spec- 
trum of the Luttinger chain does not contain low-lying 
singlets with fc = tt and even parity under spatial- 
reflections, yielding to a critical exponent X > 1 for the 
operator ||2Jl, and therefore proving its irrelevance. This 
is also illustrated in Fig. ^ where we report exact di- 
agonalization results in the weak frustration regime of 
the isotropic chain^ However, we stress that our con- 
clusions follow directly from the bosonization technique 
and represent an exact result for the isotropic spin-half 
Heisenberg chain. 

The irrelevance of the n.n.n. spin-Peierls perturbation 
in the regime of weak frustration can be also conflrmed 
numerically from the Lanczos calculation of the associ- 
ated generalized susceptibility. As explained in Ref . 
this can be done by adding to the Hamiltonian the con- 
sidered perturbation rjO and then computing the ground- 
state energy per site e{rj) for few values of 77 in order to 
estimate numerically x = —d^^{v)/d'r]^\rj=o as the limit 
X = Hm^^o 2(eo — e{'q))/rf. For instance, for the 28-site 
chain at J2/J1 = 0.1, this procedure gives Xnn — 6.4 
and Xnnn — 10^^, for the n.n. and n.n.n spin-Peierls 
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FIG. 3: Upper panel: singlet (open dots) and triplet (full 
dots) gaps versus J2/J1 for 5 = 1.0, L = 12 (triangles), 18 
(squares), and 24 (circles). The arrow indicates the extrapo- 
lated value of J2 / Ji at the level crossing. Lines are guides for 
the eye. Lower panel: size scaling of J2 / Ji at the level cross- 
ing for (from below) 5 = 0, 0.2, 0.4, 0.6, 0.8, 0.9, 1.0. Lines are 
quadratic fits. 

operator, respectively. This fact indicates that the re- 
sponse of the system to a n.n.n. dimerization is more 
than six orders of magnitude smaller than the response 
to the standard n.n. dimerization. 

Entering the regime of strong frustration, the effect 
of the n.n.n. alternation is to reduce the spin gap A 
as it is shown in Fig. |21 at J2/J1 = 0.5. As it has 
been already pointed out^ii^ at the Majumdar-Gosh point 
(J2/J1 = 0.5) the presence of the n.n.n. spin-Peierls 
dimerization does not alter the ground-state manifold 
which is spanned by the two well-known dimerized states, 
with energy Eq = — 3/8A^, independent of S. Neverthe- 
less, the excited states are not known exactly so that 
we have calculated the spin gap using the Lanczos tech- 
nique. By extrapolating the finite-size data according 
to the expected law^ii A(L) = Aqo + a/L'^ + b/L"^, we 
have obtained the S dependence of the thermodynamic 
gap, which is shown in the same figure. Our results 
are in good agreement with the variational calculation of 
Refs. Illlalia and indicate a 8% spin-gap reduction from 
the Majumdar-Gosh model for (5 = (A = 0.235) to the 
so-called saw-tooth chain for S = 1 {A — 0.216)ii^ 

The reduction of the spin gap in the dimerized phase 
induced by the n.n.n. bond alternation leads to an 
increase of the critical value {J2/Ji)c for the fiuid- 
dimer transition. As originally suggested by Castilla 
and coUaboratorsr^the phase boundary between the two 
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FIG. 4: Critical J2/J1 ratio for the fluid-dimer transition 
versus b. 

phases can be calculated very accurately by exploiting 
the change in the excitation spectrum occurring at the 
critical point. In fact, below the critical point the system 
displays quasi-long-range antiferromagnetic order: in this 
case, gapless excitations are obtained by creating low- 
energy spinous both in triplet (spatially even) and sin- 
glet (spatially odd) states, the former having the lowest 
energy on finite size. In the dimerized phase, instead, 
the triplet excitations acquire a gap so that the ground 
state remains two-fold degenerate in the singlet sector, 
leading to the dimerized ground state which breaks the 
translational invariance. As a result, the critical coupling 
{J2/ Ji)c can be determined by performing a size scaling 
of the value J2/ Ji where the level crossing between the 



low-lying singlet and triplet occurs (see Fig.O. With this 
technique, using exact diagonalizations up to 30 sites, 
it is easy to obtain very accurate determinations of the 
critical value {J2/ Ji)c- For instance, in absence of n.n.n. 
bond alternation, we obtain a critical value {J2/Ji)c hi 
agreement with previously given estimates^iiS The re- 
sulting phase boundary for the fluid-dimer transition is 
shown in Fig. 01 indicating that for 5 — 1 the critical 
value of the frustration ratio for the fluid-dimer transi- 
tion is only increased by about 1%. 

In conclusion, we have studied the effects of a next- 
nearest-neighbor spin-Peierls dimerization on the ground 
state of the spin- half J1 — J2 Heisenberg chain, using con- 
formal field theory and exact diagonalization. We have 
found that, despite the conclusions of a recent study 
based on bosonization and renormalization group ap- 
proachesji the bond alternation represents an irrelevant 
perturbation in both the regime of strong and weak frus- 
tration. In particular, as predicted by Sarkar and Sen;^ 
the perturbation does not change the nature of the exci- 
tation spectrum, leading only to a weak reduction of the 
triplet gap in the dimerized phase, and, therefore, to a 
slight enhancement of the critical value {J2/ Ji)c associ- 
ated to the fluid-dimer transition. 
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